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ABSTRACT 

We analyse SAURON kinematic maps of the inner kpc of the early-type (Sa) barred 
spiral galaxy NGC 5448. The observed morphology and kinematics of the emission- 
line gas is patchy and perturbed, indicating clear departures from circular motion. 
The kinematics of the stars is more regular, and display a small inner disk-like system 
embedded in a large-scale rotating structure. We focus on the [O III] gas, and use a 
harmonic decomposition formalism to analyse the gas velocity field. The higher-order 
harmonic terms and the main kinematic features of the observed data are consistent 
with an analytically constructed simple bar model. The bar model is derived using 
linear theory, considering an m = 2 perturbation mode, and with bar parameters 
which are consistent with the large-scale bar detected via imaging. We also study 
optical and near infra-red images to reveal the asymmetric extinction in NGC 5448, 
and we recognise that some of the deviations between the data and the analytical 
bar model may be due to these complex dust features. Our study illustrates how 
the harmonic decomposition formalism can be used as a powerful tool to quantify 
non-circular motions in observed gas velocity fields. 

Key words: galaxies: bulges - galaxies: spiral - galaxies: evolution - galaxies: for- 
mation - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

Dynamical studies of spiral galaxies often make the distinc- 
tion between the bulge and the disk: the bulge is associated 
with the hot "spheroidal" component which mainly contains 
stars, and the disk is a rapidly rotating thin structure which 
also contains a substantial fraction of the total amount of 
interstellar gas. Line-of-sight velocity distributions are effi- 
cient probes of the dynamical structures of these systems, 
and can be used to derive the mass distribution, structural 
properties, and perturbations of the gravitational potential 
linked to m — 1 or m = 2 modes, or to external triggers 
such as interactions. These can be obtained with classical 
long-slit spectrography, although it requires a-priori assump- 
tions for the orientation of the slits. A time-consuming way 
to obtain velocity fields through single slits is to spatially 
scan the galaxy with successive individual exposures (e.g., 
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Statler 1994; Ohtani 1995). Integral field spectrography in 
the optical (IFS, Adam et al. 1989; Afanasiev et al. 1990; 
Bacon et al. 1995; Davies et al. 1997; Garcfa-Lorenzo, Ar- 
ribas & Mediavilla 2000; Roth, Laux & Heilemann 2000; 
Bacon et al. 2001; Bershady et al. 2004) was designed to 
uniformly cover the field of view which makes it possible 
to study these systems in a much improved way: the anal- 
ysis is not constrained by a priori choices of slit position, 
and the gaseous and stellar components can be observed 
simultaneously. Optical IFS is often used complementarily 
to Hi and Fabry-Perot scanning interferometers, (e.g., van 
Gorkom et al. 1986; Plana & Boulesteix 1996) as the latter 
usually have excellent spatial coverage but a limited spectral 
domain, only probing a couple of lines (mostly in emission). 

Although circular rotation is the dominant kinematic 
feature of the disk component, previous observational as 
well as theoretical studies have shown the presence of non- 
circular motions (e.g., Freeman 1965; Combes & Gerin 1985; 
Shlosman, Frank & Begelman 1989; Athanassoula 1994; 
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Shlosman & Noguchi 1993; Moiseev 2000; Regan & Teuben 
2004; Wong, Blitz & Bosnia 2004). Cold disk-like systems are 
more prone than bulges to perturbations and instabilities, 
the dissipative gaseous component being more responsive 
than the stellar component (e.g., Thielheim & Wolff 1982). 
Perturbations of the gas dynamics can lead to significant 
galaxy evolution via, e.g., redistribution of the angular mo- 
mentum, triggering of star formation, or building of a central 
mass concentration. 

Bars are potential actors in this context as they can 
drive mass inwards or outwards, and may participate in the 
complex process which eventually leads to the fuelling of an 
active nucleus (e.g., Simkin, Su & Schwarz 1980, Sakamoto, 
Baker & Scoville 2000). The presence of large-scale bars 
has been correlated with starbursts (e.g., Martinet & Friedli 
1997) and nuclear rings (Shlosman 1999), but only weakly if 
at all with nuclear activity (e.g., Knapen et al. 2000). This 
may not be so surprising as the fuelling of an active galac- 
tic nucleus (AGN), i.e., gas accretion onto a massive black 
hole, involves rather small spatial and short time-scales, so 
might not be related to the kpc scale bar. Inner small bars 
(Shlosman et al. 1989) may sometimes help to link the dif- 
ferent spatial scales (Emsellem, Goudfrooij & Ferruit 2003) , 
but they probably cannot serve as a universal mechanism. 
AGN activity is short-lived and occurs in the central astro- 
nomical units of galaxies, whereas current studies have fo- 
cused on long-lived phenomena influencing the central kpc 
scale. In order to properly study the influence of a bar on 
the redistribution of mass in the central kpc of a galaxy, 
we need to obtain constraints both on the source of the un- 
derlying gravitational potential and on some tracer of the 
on-going perturbations. This can be achieved by simultane- 
ously studying the stellar and gas dynamics. 

We have observed 24 early-type spiral (Sa) galaxies with 
the SAURON Integral Field Spectrograph (Bacon et al. 2001), 
mounted at the 4.2m William Herschel Telescope of the Ob- 
servatorio del Roque de los Muchachos at La Palma (de 
Zeeuw et al. 2002) . Most of these galaxies are found to have 
gas velocity fields that strongly deviate from that of a simple 
rotating disk. For the present study, we selected NGC 5448 
out of the Sa SAURON sample which shows a clear sign of the 
presence of a bar in its photometry and which has a signifi- 
cant amount of ionised-gas (Falcon-Barroso et al. 2005). In 
this paper, we analyse both the stellar and gas kinematics 
obtained with SAURON, and quantify the non-circular gaseous 
motions in NGC 5448. To this aim, we build models with 
elliptical streaming motion and compare the resulting veloc- 
ity structures with our observational data, using a harmonic 
decomposition technique. In section |21 we present our data. 
Section outlines our analysis method of the gas kinemat- 
ics. In section |1] we build a bar model for NGC 5448, and 
we present and discuss the corresponding results in sections 
Hand El 



2 THE DATA 

NGC 5448 is a barred Sa active galaxy with prominent ir- 
regular dust lanes at different spatial scales. The RC3 clas- 
sification of NGC 5448 is (R)SAB(r)a, with inclination of 
64°, systemic velocity of 1971 km s~ , position angle (PA) 
of 115°, and total B-band magnitude of 11.93. Photometry 



Table 1. SAURON data characteristics. 



Field of View 


41" X 33" 


Pixel Size 


0"8 


Instrumental Disp. 


108 km s" 1 


Spectral Range 


[4820 - 5280] A 


Spectral Features 


H/3, [Oin], Fe5015, Mg6, [Ni], Fe5270 



shows that this galaxy hosts a large-scale bar with two well- 
defined spiral arms emerging near the ends of the bar (Es- 
kridge et al. 2002). The inner parts of the arms are somewhat 
patchy and form a broken ring, the outer parts of the arms 
being smoother (Sandage & Bedke 1994). This galaxy has 
a nuclear elongated feature of about 10" (Fathi & Peletier 
2003) with bluer V — H colour than its surroundings (Car- 
ollo et al. 2002). Ho, Filippenko, & Sargent (1995) classified 
the nucleus as an AGN: its central aperture spectrum ex- 
hibits a clear broad-line emitting region. Laine et al. (2002) 
classified NGC 5448 as a non-Seyfert. 

2.1 SAURON Observations 

We observed NGC 5448 with SAURON on April the fourteenth 
2004. Detailed specifications for the instrument, data, re- 
duction procedure, and the data preparation procedure can 
be found in Bacon et al. (2001), de Zeeuw et al. (2002), 
and Emsellem et al. (2004). A brief summary of the instru- 
ment and data characteristics is presented in Table We 
obtained 4 exposures of 1800 s each, producing a total of 
1431 simultaneous galaxy spectra per frame, together with 
146 sky spectra l!9 away from the main field. We also ob- 
served standard stars to be used for accurate velocity, flux, 
and line-strength calibration of our observed galaxy. Wave- 
length calibration was done using arc lamp exposures taken 
before and after each target exposure. The standard reduc- 
tion was performed using the XSauron software package de- 
veloped at the CRAL (Lyon), providing a fully calibrated 
datacube for each individual exposure. We merged the 4 
spatially dithered exposures after a resampling to a com- 
mon spatial scale of / .'8 x 0'.'8 per pixel, leading to a total 
of 1973 spectra within our field. 



2.2 Stellar and Gaseous Kinematic Maps 

The signal-to-noise (S/N) ratio of individual spectra varies 
considerably throughout the field: in order to homogenise 
the S/N we applied the Voronoi 2D-binning technique of 
Cappellari & Copin (2003) . We have set a target S/N value 
of 60 to be able to derive reliable higher-order moments 
of the stellar velocity distribution. Penalised Pixel Fitting 
(pPXF, Cappellari & Emsellem 2004) was used to derive the 
stellar line-of-sight velocity distribution parametrised by a 
Gauss-Hermite series up to and including the fourth order 
(van der Marel & Franx 1993; Gerhard 1993). The opti- 
mal stellar template in pPXF was built with a combination 
of stellar spectra from Jones (1997) and synthesis spectra 
from Vazdekis (1999) as in Emsellem et al. (2004). The de- 
rived stellar mean velocity V, mean velocity dispersion a, 
and higher order velocity moments /13 and /14 maps, are 
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Figure 1. Top Left: Digitised Sky Survey image with SAURQN footprint and north-east orientation arrow. All other panels show the 
SAURON data. The stellar flux map and unsharp-masked SAURON image are given in mag arcsec - 2 with arbitrary zero point, and north-east 
direction as indicated. The titles are indicated at the bottom right corner of each panel with plotting ranges according to the top color 
bar. All SAURON maps have the same orientation with overplotted stellar contours in magnitude steps of 0.25, and all velocities and 
velocity dispersions are given in km s — . 



presented in bottom row of Fig. We discuss the maps in 
section |S] 

The method which we have applied to derive the 
emission-line kinematics is described and validated in Sarzi 
et al. (2005) and Falcon-Barroso et al. (2005). It consists of 
iteratively searching for the emission-line velocities and ve- 
locity dispersions, while solving at each step for both their 
amplitudes and the optimal combination of the stellar tem- 
plates over the full SAURON wavelength range. No masking of 
the regions affected by emission is thus required. 

Our binning scheme was optimised for the derivation 
of reliable stellar kinematics, but is not necessarily well 
adapted to analyse the emission-line gas. Since the gaseous 
component commonly rotates faster than the stellar compo- 
nent, to obtain reliable gas kinematics, in general a lower 
S/N is required. Moreover, gas and stars do not necessar- 
ily share the same spatial distributions. As a consequence, 
the stellar bins are often larger than necessary for the gas, 
implying a loss of spatial information for analysing the 
gaseous component (Fathi 2004). Accordingly, the data that 



we present in this paper has a minimum stellar S/N of 60 
and a minimum emission-line amplitude-over-noise of 3. 

Finally, we have carefully inspected the emission-line 
profiles to search for possible asymmetric profiles but have 
not found any significant deviation from a pure Gaussian, 
nor any significant indication for more than one component 
being present in the profiles, at least at the instrumental 
resolution of SAURON. We know from other SAURON studies of 
e.g., NGC 1068 (Emsellem et al. 2005) that we can resolve 
several components in very strong AGNs. 



3 ANALYSING GAS VELOCITY FIELDS 

3.1 Circular and Non-Circular Kinematics 

In this section we start with a simple mathematical model, 
and then add complexity as required by the observational 
constraints. In the simplest model, the galactic disk is purely 
rotational, has negligible velocity dispersion, and is infinites- 
imally thin. In the presence of axisymmetric radial and ver- 
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tical velocities, and when taking into account the effect of 
projection and the convention that positive line-of-sight ve- 
locities correspond to recession, the line-of-sight velocities 
Vies can be represented by 



Vl os (R, 1p, l) = Vsys + Vrot (R) COS l/j Sin % 

+V ra d(-R) sin ^ sini 
+V Z (R) cosi, 



(1) 



where V syB is the systemic velocity of the galaxy, V Iot and 
Kad are the rotational and radial velocities, and V z is the 
vertical velocity component, which we set to zero throughout 
this paper. The inclination i of the galaxy ranges from i = 0° 
for a face-on viewing and i = 90° for an edge-on viewing. 
Furthermore, (R, tp) are polar coordinates in the plane of the 
galaxy related to observable Cartesian coordinates x and y 
(in the plane of the sky) by 



cos ip 
sini/) ■ 



-(x - X ccn ) sin (f> + (y - Y ccn ) cos (f> 
R 

-(x - X ccn ) cos (j>o - (y - Y ccn ) sin <p 
7? cos i 



(2) 



where X ccn and Y ccn are the coordinates for the centre, and 
<po is the PA of the projected major axis of the disk, mea- 
sured with respect to north in counterclockwise direction. 
This simple model cannot explain most observed velocity 
fields. The gas kinematics in real galaxies exhibits radial 
and/or vertical motions due to, e.g., the presence of bars, 
spiral arms, which create angle-dependent velocities which 
cannot be explained using Eq. Q, and hence additional in- 
gredients in the analysis method are required. 

Several attempts have been made to investigate more 
complex velocity fields (Sakhibov & Smirnov 1989; Canzian 
1993; Schoenmakers, Franx & de Zeeuw 1997 [hereafter 
SFdZ]; Fridman & Khoruzhii 2003; Wong et al. 2004; Kra- 
jnovic et al. 2005). The main idea of these techniques is to 
deal separately with the "unperturbed" underlying compo- 
nent and the residual velocity field, obtained by subtracting 
the reconstructed unperturbed velocity field from the data. 
To solve the problem, as stated, one would have to know the 
unperturbed velocity component at a given radius. Unfortu- 
nately we do not know beforehand the unperturbed velocity 
field. A powerful technique for unveiling the perturbations 
makes use of the expansion of the velocity information in 
Fourier harmonic components. Following the terminology of 
SFdZ (see their appendix for the detailed derivation), given 
that the line-of-sight velocity can be expressed as a two- 
dimensional analytic function of galactocentric coordinates, 
it can be written as a Fourier series: 



1 io 



cos nip + s n (R) sin nip) sini, (3) 



where k is the number of harmonics, and c n and s n give us 
information about the nature of the perturbations, and they 
are tightly connected to the underlying potential. Note that 
ci = Vrot and si = Kad - Furthermore, numerical simulations 
by Burlak et al. (2000) have shown that this formalism is 
very stable and its results are rather insensitive to the pres- 
ence of holes in the velocity field. 



3.2 Tilted-Ring and Harmonic Decomposition 

Our quantitative analysis uses the formalism of Eq. 
combined with the tilted-ring method (Rogstad, Lockart & 
Wright 1974; Rogstad, Wright & Lockart 1976), inspired 
from the ROTCUR routine in the GIPSY package (Begeman 
1987; van der Hulst et al. 1992). Accordingly, we divide 
the galaxy disk into concentric ellipses within which we fit 
the rotation velocity and the set of geometric parameters of 
Eq. 0. The method assumes that each pixel in the velocity 
field is identified with a unique position in the galaxy and 
that circular rotation is the dominant feature. We assume 
that our measurements refer to positions on a single inclined 
disk, i.e., that we do not expect a strong warp and/or over- 
lapping spiral arms within the observed field. Accounting 
for harmonic perturbations up to order 3, this method re- 
sults in a large parameter-space problem, and it requires 
some additional assumptions to obtain a physically mean- 
ingful solution. We obtain the parameters by proceeding as 
follows: 

(i) We start from a simple inclined two-dimensional disk. 
In case of erroneous inclination, one would detect system- 
atics in the residual field as explained in section 13.31 Us- 
ing a two-dimensional disk produces a robust model, in the 
sense that warps and other deviations from this basic model 
should be prominent in the residual maps. 

(ii) Secondly, we fix the dynamical centre of the galaxy 
to coincide with the photometric centre. The photometric 
center is obtained by fixing the center of the H band image, 
since the //-band image is less affected by dust. The HST 
V-band image and the SAUR0N image are then aligned with 
the HST //-band image. 

(iii) At every radius, we fit yi os = Co + ci cos tp sin i, vary- 
ing Co, PA, and ci. The result is used to fix the systemic 
velocity (V sys = Co) simply as the mean value over radius, 
using Tukey's bi-weight mean formalism (Mosteller & Tukey 
1977), which is particularly advantageous for being non- 
sensitive to outliers. After Vsys, we fit the PA and ci, fol- 
lowed by fixing the PA in the same way as the V sys was 
fixed. 

(iv) For each ring, we import the values derived accord- 
ing to the recipe above, and fit the rotation curve and the 
higher harmonic components (up to order 3) applying a x 2 - 
minimisation scheme. These results are presented in Fig. [5] 

Analysing the first, second, and third harmonic compo- 
nents only, is an efficient way to analyse specific elements of 
the perturbations on the underlying potential. A perturba- 
tion of order m creates m—1 and m+1 line-of-sight velocity 
terms (see e.g., Canzian 1993 and SFdZ). As a result, our 
fitted ci, si, C3, and S3 terms contain information about 
possible m = 2 perturbations. 

When sectioning the field into concentric rings, the ra- 
dial extent of each ring is a free parameter which has to be 
pre-determined. For a range of possible ring radii, we have 
examined the fitted velocity information. Determination of 
the radial thickness of the rings becomes a trade-off between 
smoothness of the fitted velocity field, and robustness of the 
fit. Radially thin rings result in smooth fits, and thus smooth 
residual (Data - Fit) fields, whereas radially thick rings re- 
sult in lower errors for the fitted parameters. Patchy gas dis- 
tribution causes non-uniform distribution of points, and it 
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is important to make sure that all the rings include enough 
points to fit harmonic parameters to reasonable accuracy. 
This is particularly important in the innermost and outer- 
most points. We adopt a geometric increase of the ring radii 
(meaning that the radial width of the rings is increased by 
a factor 1 + step). Here, we adopt a step of 0.2, and only fit 
the inner 30" of the SAURON velocity field, since rings larger 
than this radius are not sufficiently covered and hence do not 
deliver reliably derived parameters. Finally, since we are si- 
multaneously fitting 6 harmonic components, we make sure 
that each ellipse contains at least 25 data points, in order 
to obtain reliable errors. 



3.3 Errors 

Warner et al. (1973) and van der Kruit & Allen (1978) 
showed that assuming wrong input disk parameters could 
cause recognisable signatures in the residual velocity fields. 
SFdZ quantified these signatures in terms of harmonic ex- 
pansion, and showed that the Fourier components due to 
non-circular motions will mix with those due to erroneous 
disk parameters X ccn ,Y ccn ,i, and PA. Errors in the kine- 
matic centre translate to the Co, C2 and S2 terms, while PA 
errors appear in the si and S3 terms, and inclination errors 
in the ci and C3 terms. In the case of erroneous kinematic 
centre, the Co and the C2 terms fall off as whereas the 
presence of a radius-independent m = 1 mode shows the 
relation Co ~ 3c2 (Schoenmakers 1999). Our chosen centre 
coordinates do not deliver any of these situations. Another 
source of error is "pixel sampling" for which errors propa- 
gate onto the ci (as the rotation curve rises more gradually) 
and the C3 term (as if inclination changes). This effect may 
also propagate onto the si and S3 parameters, although not 
significantly (Wong 2000). Our data and analysis does not 
require consideration of this effect since our spatial resolu- 
tion is high enough. 

Our data extraction code delivers error measurements 
for each individual pixel (Falcon-Barroso et al. 2005). We 
calculate the errors for the tilted-ring and harmonic param- 
eters by means of Monte Carlo simulations. Repeated appli- 
cation of the tilted-ring method to the Gaussian randomised 
gas velocity field yields the uncertainties on the harmonic pa- 
rameters. Our 500 simulations show that choosing geomet- 
rically increasing ring radii with steps of 0.2 indeed yields 
satisfactory errors. The errors are low, and the parameter 
profiles are represented by an adequate number of points. 



4 A BAR MODEL FOR NGC 5448 

The SFdZ formalism is a very powerful tool for analysing 
velocity fields. One can analytically derive the higher-order 
harmonic parameters for potentials of desired form. Schoen- 
makers (1999) derived the harmonic parameters for a lop- 
sided potential, and a more extensive expansion was done by 
Wong et al. (2004), who derived the higher-order harmonics 
for simple bar and two-armed spiral perturbations. 

We import their model for a globally elongated grav- 
itational potential, which, in the plane of the galaxy, con- 
sists of an axisymmetric component $o(-R) from Kuijken & 
Tremaine (1994), with a bi-symmetric potential of the form 



$2 (-R) cos 2$ added to it. In the case of flat rotation curve, 
the potential is written as 

$o(-R) = v 2 c lnfl, 

(4) 

$2{R) = -evl I 2, 

where v c is the circular velocity and e is the elongation of 
the potential. Introducing a constant damping term into the 
equations of motion and considering only the m — 2 pertur- 
bation term (Lindblad & Lindblad 1994; Wada 1994), linear 
epicyclic approximation of Lin & Shu (1964) was applied to 
obtain an analytical solution for gaseous orbits in this barred 
potential. Since this model is restricted to m = 2 perturba- 
tion, it affects only the first and third harmonic parameters 
(ci, Si, C3, and S3). These depend on the amplitude of the 
damping term (A), the ellipticity of the potential (e), the 
corotation radius (CR), and the viewing angle (9), where 
— corresponds to end-on view. 

We use these analytic harmonic parameters to recon- 
struct the signatures of prominent non-circular motions in 
the observed velocity fields. We build a library of models 
with varying bar parameters (0 < e < 0.5, < A < 0.5 
and a range of bar sizes), viewed from different angles 
(0° < 9 < 180°). The effects of A and e on the velocity 
field are predictable. The bar signature is weaker for a larger 
damping factor and/or a smaller flattening of the potential. 
However, the effect of varying 9 is not straightforward. 

Investigating our library of dissipative bar models, we 
find that only for a certain combination of parameters the 
model exhibits a strong twist of the gas zero- velocity curve 
similar to our data (Fig. 0. We quantify the comparison 
by calculating, for each bar model in our library, the dis- 
crepancy between the reconstructed velocity field and the 
observed velocity field, in terms of the goodness-of-fit pa- 
rameter x 2 - Care has to be taken at the inner Lindblad res- 
onance (ILR), since there the linear epicyclic theory breaks 
down. At the ILR, due to the presence of non-linear terms, 
which are not accounted for by the analytic epicyclic the- 
ory, the weak-bar model is not expected to reproduce the 
velocity field. Hence, for a small region around the ILR, we 
interpolate the bar models. 

The minimum goodness-of-fit parameter Xmm yields the 
best-fit bar model, for which we show in Fig. |2] the recon- 
structed velocity field and harmonic terms (red filled circles) . 
The best-fit model accounts for an m = 2 perturbation of a 
single bar, with potential ellipticity e — 0.15 ± 0.02, damp- 
ing term A = 0.12 ± 0.03, viewing angle 9 = 19 ± 3°, and 
CR = 37 ± 4". The 3a errors on the bar parameters follow 
from the bar models for which the difference between the 
corresponding \ 2 anc ' Xmin is below the 99.7% level. From 
the harmonic terms of these bar models, we also obtain an 
estimate on the uncertainty in ci, s\, C3 and S3, indicated 
by the red error bars in Fig. [5] Given the significant sec- 
ond harmonic term and the simplicity of the analytic bar 
model, it is not surprising that it cannot provide a perfect 
fit to the data. Still, this generic bar model does reproduce 
the main features in the observed velocity field, as well as 
the overall behaviour of the first and third harmonic terms. 
This supports the case of a bar as the main driver behind 
the observed non-circular motions. 
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Figure 2. Top row, from left to right: observed SAURQN gas velocity field of NGC 5448; reconstruction based on the harmonic decomposition 
of the SAURQN gas velocity field; residual field (data - harmonic reconstruction); and the analytic bar model, which reproduces best the 
main kinematic features of the observed gas velocity field. Middle and bottom row: the harmonic parameters as a function of outer radius 
of each ring, and v* = cisini. The over-plotted red filled circles are the analytically calculated first and third harmonic terms for the 
bar model (the second terms are zero by construction), with the error bars corresponding to the 3c (99.7%) confidence level (see section 
4 for details). The orientation of the maps is the same as in Fig.Q and the mean PA value is given in the north-east direction. 



5 RESULTS 

The gas velocity field decomposition of section 13.21 pro- 
vides the radial profiles for the primary disk parameters and 
the higher order harmonic terms. Assuming the inclination 
i — 64° from the RC3 catalogue (de Vaucouleurs et al. 1991), 
we apply this method to our observed gas velocity field in an 
iterative manner by first varying the V S ys and PA. We then 
fix these parameter to their mean values and we proceed by 
deriving the V ro t and higher harmonic terms (see Fig.[2j. We 
detect a varying zeroth harmonic term (V sys ) with variation 
amplitude up to 40 km s _1 in the central 10". We find a 
strong PA twist of about 30° in the central 10", and out- 
side this radius the PA goes back to the central values. The 
circular velocity component ci rises steeply and peaks at 
5" radius, followed by a very slow decline out to the outer 
radii. We find that the behaviour of the si term could be 
compared with the Vsys or the PA: the si term changes sign 
at around 12", and at the same galactocentric radius, the 
V sys and the PA change sign around their mean values. The 
si exhibits amplitudes of order of 0.2cisini. The bottom 
row in Fig.|2|shows that, apart from the S2 term, the higher 
harmonic terms significantly deviate from zero, and indeed 



account for the considerable deviation of the velocity field 
of NGC 5448 from pure rotation. 

In the following sections, we describe the distinct fea- 
tures that we find in our kinematic maps. 

5.1 The Stellar Component 

The SAURDN flux map in Fig.^displays a smooth stellar dis- 
tribution and the presence of prominent dust lanes to the 
south of the nucleus. The stellar kinematics shows a promi- 
nent regular disk-like rotation with a smaller inner stellar 
disk within the central 7" radius. The stellar velocity disper- 
sion decreases towards the centre, and at the location of the 
maximum line-of-sight velocities for the central structure, 
hs anti-correlates with V which supports the argument that 
the central component of NGC 5448 is a central disk. We 
approximate the stellar velocity field with that of an expo- 
nential thin disk (Freeman 1970) to emphasize the kinematic 
signatures of the central disk. Fig. [3] shows this simple in- 
clined disk model, where the inner stellar disk rotates faster 
than the outer disk. The best-fit model for the outer stellar 
disk, omitting the central 7" and assuming an inclination 
fixed to the RC3 catalogue value, yields V / sys = 2002 km s -1 , 
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a disk scale length Rd = 18, and a disk PA= 91°. Repeating 
the same exercise for the region inside 7", we find that the 
central disk is misaligned with respect to the outer disk by 
~ 13°. 



5.2 Gas Distribution and Kinematics 

Fig-Qdemonstrates that along the strong dust lanes, the gas 
shows a patchy distribution, with an asymmetric elongation 
of [O ill] gas towards the east as well as the galactic poles. 
The H/3 distribution is found to be more regular, and we 
have found that all kinematic features of the [O ill] gas are 
accompanied by H/3 emitting gas. The [O ill] /H/3 map in 
Fig. displays a prominent ring-like structure at 5"- 15", 
indicating high ionisation of the gas in this region. 

The gas velocity map clearly shows very prominent 'S'- 
shaped zero-velocity curve with very sharp edges indicating 
very strong non-circular gas motions (Peterson & Huntley 
1980). This 'S '-shaped gaseous zero- velocity curve is con- 
firmed by the PA variation derived by the tilted-ring decom- 
position in Fig. |21 Zero-velocity curve twists can be indica- 
tive of strong radial motions, created by, e.g., interactions, 
mergers or elliptical streaming due to a barred potential. 
Wong (2000) distinguished differences between radial mo- 
tion mechanisms by comparing the S3 versus si harmonic 
terms, and found that the behaviour of S3 versus si is differ- 
ent between a bar model and a pure axisymmetric or warped 
model. As demonstrated in Wong et al. (2004; Fig. 5) the S3 
versus si of a warped disk lie on a positive slope. Moreover 
in case of an externally triggered radial flow, the points in 
the S3 versus Si graph should lie on a zero-slope. Setting up 
the same diagnostics, in Fig. 01 we find that for our data, 
the S3 versus si curve has a negative slope very similar to 
that of the bar model presented in Fig. [5] This confirms the 
signature of elliptical streaming in the gas velocity field of 
NGC 5448. 

Comparing the stellar maps with the gas maps, we find 
that the gas velocities are considerably higher: This is con- 
firmed with the rotation curve presented in Fig. |H] We find 
that the gaseous a map displays some features comparable 
to the stellar a map. Although the stellar velocity dispersion 
is more steeply rising toward the centre, both maps show a 
prominent dip in the centre. This supports the hypothesis 
of the presence of a central dynamically cold disk-like struc- 
ture, present both in the stellar and gaseous component. 
Both maps also show a dispersion increase in bi-polar di- 
rections just a few arcseconds from the centre. This is much 
more prominent in the gas map, as its a reaches values ~ 250 
km s _1 . This may be indicative of significant outflows in this 
AGN host. 



5.3 Position Angle Values 

The projected and deprojected PA profiles from Laine et 
al. (2002) in Fig. display strong variations in the central 
few arcseconds. The HST images display asymmetrically dis- 
tributed dust within the central few arcsecoonds, which may 
cause these strong central PA variations (also found by Ko- 
rnreich et al. 2001). At larger radii, we find that the PA 
profile shows a variation of the order of 10°, but overall 
it stays fairly constant around the nominal RC3 value, i.e. 
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Figure 4. First and third sinusoidal harmonic terms normalised 
to v* = ci sini. Both the data and the bar model show a negative 
slope (dashed and solid line respectively), indicating that elliptical 
streaming, due to the bar potential, is the main driver of the radial 
motions observed in the gas velocity field of NGC 5448. 



115°. Our tilted-ring decomposition of the gaseous velocity 
field, yields a mean gas kinematic PA of 119° ± 5, which is 
consistent with the photometric PA from RC3. In the case 
of the large stellar disk, as described in section l5~T1 we find 
an offset of ~ 25° between the stellar kinematic PA and the 
photometric mean PAs. This misalignment could be due to 
the effect of the bar on the observed stellar velocity field. 



5.4 The Inner Few Arcseconds 

In Fig.HJwe zoom into the central few arcseconds of the stel- 
lar and the gaseous velocity fields of NGC 5448 to emphasize 
the observed kinematic differences between the two compo- 
nents. The gaseous component exhibits a strongly twisting 
zero- velocity curve, and the gas rotates faster than the stars. 
We find that the inner stellar disk not only is misaligned with 
the outer stellar disk, but also is shifted to the south-east. 
In contrast, we find from the photometric profiles in Fig. [H] 
that the ellipticity of the central disk is similar to that of 
the outer disk (~ 0.55), as well as their PAs. The promi- 
nent dust lanes at the south of this galaxy may explain the 
kinematic misalignment between the outer and inner stellar 
disks. 

In the central ~ 4" the gas isovelocity contours are com- 
paratively straight and parallel. In the same region, the stel- 
lar isovelocity contours are indicative of the presence of a 
central disk. The tilted-ring decomposition of the gas veloc- 
ity field provides the gas rotation curve illustrated in Fig. |S| 
The gas rotates about 70 km s _1 faster than the stars in 
the region outside ~ 4", with some indication that outside 
this radius the negative velocities reach lower values than 
the positive line-of-sight velocities. This asymmetry could 
be caused by complex dust distribution or attributed to a 
signature of lopsidedness. 
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Figure 3. A thin isothermal disk model for the stellar velocity field of NGC 5448. The circle marks the 7" region within which we find 
a disk-like structure. The disk model, fitted to the field outside this region, implies for the outer disk a scale length of 18", V S ys= 2002 
km s~ , and PA= 91°. The orientation of the maps is the same as in Fig. HI 
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Figure 6. Zooming into the central few arcseconds of the stellar and gaseous velocity maps of NGC 5448, using the same velocity range 
as in Fig. Indicated are the north-east direction (arrow), the photometric PA (straight line) and the photometric centre (cross). The 
over-plotted circle indicates the 7" radius for comparison with Fig.|3] In the right panel, we present the stellar rotation curve (extracted 
along the photometric PA) together with the gas rotation curve derived from tiltcd-ring decomposition. 



5.5 Dust Effects or Lopsidedness? 

Although the unsharp-masked SAURON image shows strong 
signature of dust-lanes in the region south of the nucleus, 
the central dust distribution is easier to study when looking 
at the HST images and the V — H colour map in Fig. |S] 
The colour map shows a strong dust lane about 1" at the 
south-east of the galaxy nucleus as well as an asymmetri- 
cally distributed overall dust distribution. The larger scale 
WFPC2 image shows that this central dust lane is accom- 
panied by several other dust lanes further away from the nu- 
cleus. Using the Galactic extinction law of Rieke & Lebofsky 
(1985), the HST images yield an average extinction value of 
Av — 0.5 over the central 3", with a maximum value of 
1.5. The extinction increases significantly towards the cen- 
tre, with a prominent peak a few arcseconds south of the 
nucleus (see Fig. |5J. 

Some of the observational effects of the asymmetrical 
dust distribution could be interpreted as signatures of lop- 
sidedness. Lopsidedness in galaxies has been investigated by 
e.g., Baldwin, Lynden-Bell & Sancisi (1980) and Swaters 
et al. (1999). The amplitude of this effect may depend on 



galactocentric radius, and viewed from different viewing an- 
gles will produce different signatures in the observed velocity 
field. These studies have shown that the residuals are domi- 
nated by the zeroth harmonic term for viewing angle of 90° , 
and by the second harmonic term for viewing angle of 0° . In 
the special case of radius-dependent lopsidedness Co ~ 3c2- 
NGC 5448 has a large inclination, thus the zeroth term dom- 
inates, and we find no correlation between the Co and the C2 
terms. Our harmonic decomposition results display strongly 
varying zeroth and second terms. This was shown by SFdZ 
to be a poassible indication of an m — 1 perturbation, i.e. 
lopsidedness (c.f., Fig 2. in SFdZ). In this case, subtracting 
the rotational component from the perturbed velocity field, 
one should find a ring-like feature. In Fig. |7| we present this 
residual field for NGC 5448 and find no ring-like feature. 
Therefore dust is most likely the cause for the asymmetric 
features that we observe in NGC 5448. 
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Figure 5. The top panels show HST/NICMOS //-band (F160W) 
and WFPC2 V-band (F606W) images of NGC 5448. The bottom 
left panel is the V — H map obtained by convolving the //-band 
image with the V-band image PSF and vice versa, and subtract- 
ing one from the other. All maps have pixel-size of 0'.'0455, and 
the same orientation as the DSS image with the centre marked 
by a cross. The bottom right profiles are the observed (open sym- 
bols) and deprojected (filled symbols) ellipticity and PA profiles 
from Laine et al. (2002). Circles represent values derived from the 
NICMOS image, and triangles show the values derived using the 
DSS image. Here we only focus on the region beyond 1". 
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Figure 7. Residual field after subtracting the best-fit rotational 
model from the observed gas velocity field. This map reveals the 
strong non-circular motions in NGC 5448, however, no ring-like 
feature is found. 
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Figure 8. Observed (left) and deprojected (right) DSS image of 
NGC 5448. The inclination of 64° and PA of 115° from the RC3 
catalogue have been used to deproject the image according to the 
two-dimensional deprojection of a circular thin disk. 



5.6 The Large-scale Bar 

In Fig |5J we find that the harmonic parameters derived 
from the data are consistent with the bar model outlined 
in section 2] Our simple model is limited to an m = 2 
perturbation. Complicated flow patterns in bars studied by, 
e.g., Lindblad (1999), imply a significant contribution from 
m — 4 and higher-order modes. The third harmonic terms 
are also influenced by m = 4, which could explain the dif- 
ferences between our model bar and the data for these pa- 
rameters, but a detailed model is beyond the scope of this 
paper. 

We now compare the modelled bar with the deprojected 
Digitised Sky Survey (DSS) image. The inclination from the 
RC3 and the PA from our tilted-ring model were used to 
deproject the DSS image in a simple way by assuming that 
the galaxy disk is thin. The outer spiral arms are clearly 
emphasized in the deprojected DSS image (Fig. 0. 

It is somewhat difficult to identify the size of the large- 
scale bar from the image. From our bar model we find a 
CR of about 37". In practice, this is done by assuming that 
at CR the circular frequency of the axisymmetric potential 
is identical to bar angular frequency. According to this def- 
inition and since we use a flat rotation curve, the ILR is 
located at 1 - l/y/2 of the CR, i.e., at about 11". If we as- 
sume that the bar in NGC 5448 ends close to the starting 
point of the spiral arms (Sanders & Huntley 1976), we can 
also then associate its CR with the inner radius of the arms. 
Fig. |H| shows that our analytic bar model is of reasonable 
size. The assumed CR radius is confirmed by the ellipticity 
profile presented in Fig. |S] since the high ellipticity plateau 
starts at around this radius. At around the ILR (i.e., 11"), 
there is a good agreement between the si profiles from the 
data and the model, with in both cases a change in sign at 
the ILR. 



6 DISCUSSION AND CONCLUSIONS 

We have used linear theory to construct a bar model (m = 
2 mode) for NGC 5448. The best-fit model considers the 
effects of a single weak bar with a potential ellipticity e = 
0.15 and a damping amplitude of A = 0.12. We compare the 
harmonic terms of the model with those derived from the 
data and find an overall consistency (see Fig. [Sj. Using this 
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model, we have been able to associate the radial motion of 
gas with that of the large-scale bar (see Fig. 

NGC 5448 exhibits clear signatures of the presence of 
other components than a single bar, which affect the ob- 
served velocity field. Inspecting the photometry and the cen- 
tral parts of the gaseous and stellar velocity field, we have 
detected a central rotating disk-like stellar component em- 
bedded in the larger disk (see Fig. EJ- The WFPC2 image 
shows one very sharp dust lane 1" from the centre of this 
galaxy. This image also shows three more fuzzy and almost 
parallel dust lanes further out in the south-east direction. 
The dusty centre is also apparent in the V — H image in 
Fig. Q>] Dust is inhomogeneously distributed and continues 
down to the very bright nucleus. Light from the central light 
source is absorbed by the dust asymmetrically and the cen- 
tre appears to be located at the north-west of its actual 
position. 

Investigating the projected and deprojected elliptic- 
ity profiles derived from the //-band data from Laine et 
al. (2002) confirms the presence of a central disk in the in- 
ner 7". Although the projected ellipticity profile decreases 
at ~ 4" and outwards, we cannot pin down the exact size 
of the central disk. The observed ellipticity decrease could 
be partly caused by the strongly asymmetric dust distribu- 
tion. The stellar kinematic maps show that the central disk 
rotates faster than the main disk, and our observed gas dis- 
tribution and kinematics indicate that this central disk also 
hosts gas which rotates faster than its stellar counterpart. 
This is not unexpected, since due to its disspative nature, 
the gas is not slowed down by asymmetric drift. It is impor- 
tant to note that, in NGC 5448, the gas velocities are well 
ordered but larger in magnitude than the stellar velocities. 
However, Fig. |S| indicates that this velocity difference dies 
out at larger radii. The difference between the stellar and 
gas velocities can be easily explained as due to the fact that 
the stars are in a thicker structure (bulge) in the central re- 
gions, while at larger radii both the gas and the stars are in 
a flatter disk, and both move closer to the circular velocity. 

It is known that bars are efficient in transferring mass 
towards the inner regions of their host galaxies. The cen- 
trally concentrated matter may be able to form a central 
disk (Yuan & Yen 2004), or an inner bar (e.g., Maciejewski 
& Sparke 2000, Englmaier & Shlosman 2004). Our analysis 
has shown that the non-circular gas kinematics in NGC 5448 
could be driven by the large-scale bar. The central disk could 
then have been formed as a result of the gas accumulation 
at the centre. 

To conclude, we have been able to analytically model 
the bar signatures in the AGN host NGC 5448. We have un- 
veiled a central disk, and have distinguished the effects of a 
lopsided perturbation from strong dust features. NGC 5448 
hosts considerable amounts of dust which is asymmetrically 
distributed all the way to the centre, resembling the "Evil 
Eye" galaxy (Braun, Walterbos & Kennicutt 1992). This 
study shows the power of the harmonic decomposition for- 
malism to quantify non-circular motions in observed velocity 
fields, and we plan to apply this approach to the full set of 
bulges observed with SAURON. 
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